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Abstract 

We study dynamics of a QGP fluid induced by energetic partons propagating through it. We 
construct a (3+l)-dimensional QGP-fluid+Jet model. When a jet traverses a uniform fluid, it 
induces a Mach cone structure of energy density distribution and a vortex ring surrounding a 
path of the jet. When a pair of jets travels through a radially expanding fluid, low momentum 
particles are dominantly induced at large angles from the quenched jet. This result is quahtatively 
consistent with observation of the CMS Collaboration at LHC. 



1. Introduction 

Quark gluon plasma (QGP) created in high-energy heavy-ion collisions behaves like a perfect 
fluid. Energetic partons created through initial hard scatterings are subject to lose their energy 
during traversing the QGP medium, which is called jet quenching 1 1, 2, 3]. Thus the jet quench- 
ing phenomenon is a promising tool to extract information about stopping power of the QGP. 
These partons could induce collective flow in the QGP by depositing energy and momentum. In 
fact, the CMS Collaboration observed a large number of low momentum hadrons at large angles 
from a quenched jet |4|]. The sum of these low momentum hadrons together with a quenched 
jet balances that of the momentum of a jet created at the same position and propagating in the 
opposite direction. Therefore the result can be interpreted as manifestation of a wake of the QGP 
fluid by the jet quenching. 

In this paper, we construct a relativistic hydrodynamic model to demonstrate a response of 
the QGP fluid to jet propagation through it. We model a source term in relativistic hydrody- 
namic equations which exhibits deposit of the energy and momentum from traversing jets. With- 
out linearisation we solve these non-linear hydrodynamic equations numerically in full (3-1-1)- 
dimensional space and describe responses of expanding QGP medium to propagation of energetic 
jets. 

2. QGP-Fluid+Jet Model 

To study the dynamics of the QGP fluid induced by jets, we use relativistic hydrodynamic 
equations with source terms 

5^r^^(x) = r{x). (1) 

Here T^^'' is the energy-momentum tensor of the QGP perfect fluid and is the density of the 
four momentum deposited from the traversing jets. We employ massless ideal gas equation of 
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state, P — e/3, where e and P are energy density and pressure, respectively. We assume that 
the deposited energy and momentum are instantaneously thermalised inside a fluid cell. For 
simplicity, we consider collisional energy loss only. When a massless jet partcle travels through 
the QGP fluid, the source term is given by 

J\x) = -^<J*" [x - x,,,{t)) , J{x) = ^f{x). (2) 

We solve Eq. ([l) numerically without linearisation for a given initial condition. In this way, we 
describe collective flow induced by jets on top of an expanding background QGP fluid in high- 
energy heavy-ion collisions. Particle distributions can be calculated using Cooper-Frye formula 
list] in an isochronous hypersurface in the calculation frame. 

3. Results 

We consider two different situations. As a test case, we consider one jet particle traveling 
through a uniform fluid. As a more realistic case, we consider a pair of jets traveling in opposite 
directions through an expanding fluid. We assume that jet particles are massless and travel in a 
straight hne. 

3.1. One jet traveling through a uniform fluid 

Figure [T] shows the energy density distribution and the flow velocity field of the QGP fluid 
in a plane including a jet axis at a fixed time t - 9 fm/c. Initial energy density of the fluid 
is 0.019 GeV*. The jet particle starts its propagation in the jc-direction at the speed of light 




Figure 1 : Contours of the energy density distribution and the flow velocity field of the fluid at f = 9 fm/c. 



from the origin at f = fm/c. One sees the so-called Mach cone, namely a "V-shaped" cone- 
like structure of higher energy density region in three dimensional space. The jet particle is 
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followed by a strong collective flow in the positive x direction, which leads to lower energy 
density behind the cone than that of the backgrounds. Notice here that the background energy 
density is e = 0.019 GeV"*. On the surface of the cone, flow velocity is nearly perpendicular 
to the surface. Behind the cone, one sees vortices around {x,y) - (5 fm, +2 fm) which form a 
vortex ring around the path of the jet in three dimensional space. 

3.2. A pair of jets traveling through an expanding fluid 

Next, we consider a pair of jets traveling through a QGP fluid which expands radially. Initial 
energy density distribution is 



e{t = 0, x) 



pexp(-^), 



(1^1 < R) 
{\x\ > K) 



(3) 



Here, we choose To - 0.5 GeV, cr - 1.5 fm and R - 6 fm. Without jets, this fluid expands 
isotropically. A pair of jets is supposed to be created at {x,y,z) - (3 fm, Ofm, Ofm) with the 
common initial energy 150 GeV, both of which are traveling in the opposite direction. One jet 
goes in the positive x direction and is observed as a leading jet. The other goes in the negative x 
direction, loses more energy and is observed as a sub-leading jet. 

To see the momentum acquired by the QGP fluid from the jets, we calculate total momentum 
distribution as a function of polar angle measured from the leading jet axis subtracted by the 
distribution in the case without jet propagation: 
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where 9 is the polar angle, p the momentum of the particle from the fluid and / the index of the 
momentum regions. Summation of the particles is taken over in the three different momentum 
regions such as "high-p" {p > 8 GeV), "mid-/?" (4 GeV < /:> < 8 GeV) and "low-/?" (0 GeV < 
p < A GeV). We calculate Eq. (|4]l through the Cooper-Frye formula and do not add final jets' 
momentum; 
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Fi sure l2l shows A ( fiQQ^g \ "is a function of cos 6 for three different momentum regions. Note that 
cos 6 - -I corresponds to the direction of the sub-leading jet. A single prominent peak appears 
in this direction for each momentum region and no additional peak away from the direction of 
the sub-leading jet, which is expected from the Mach-cone structure, is seen. At large angles 
from the sub-leading jet (cos 6 > -0.94), low momentum particles are dominant. This result is 
qualitatively consistent with the CMS data [4] and indicates importance of medium responses to 
jet propagation to understand the modification of jet structure. 



4. Summary 

In this paper, we studied the dynamics of the QGP fluid induced by jet propagation. We con- 
structed a QGP-fluidH-Jet model and performed two simulations. In the case of one jet traveling 
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Figure 2: The angular distribution of the fluid's momentum received from the jets. Solid, dashed and dotted lines 
correspond to three different momentum regions, "high-p" (p > 8 GeV), "mid-p" (4 GeV < p < 8 GeV) and "low-p" 
(0 GeV < p < 4 GeV), respectively. 



through a uniform fluid, a Mach cone appears in the energy density distribution and a vortex ring 
is formed around the path of the jet behind the cone. In the case of a pair of jets traveling through 
an expanding fluid, we showed that low momentum particles are dominant at large angles from 
the subleading jet. These results indicate that the QGP-fluid-i-Jet model describes the CMS data 
qualitatively and importance of the QGP medium responses to the jet propagation. 
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